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ABSTRACT 

We present the first search for spinning dust emission from a sample of 34 Galactic 
cold cores, performed using the CARMA interferometer. For each of our cores we use 
photometric data from the Herschel Space Observatory to constrain N h, T^, an d 
Go- By computing the mass of the cores and comparing it to the Bonnor-Ebert mass, 
we determined that 29 of the 34 cores are gravitationally unstable and undergoing 
collapse. In fact, we found that 6 cores are associated with at least one young stellar 
object, suggestive of their proto-stellar nature. By investigating the physical conditions 
within each core, we can shed light on the cm emission revealed (or not) by our 
CARMA observations. Indeed, we find that only 3 of our cores have any significant 
detectable cm emission. Using a spinning dust model, we predict the expected level 
of spinning dust emission in each core and find that for ah 34 cores, the predicted 
level of emission is larger than the observed cm emission constrained by the CARMA 
observations. Moreover, even in the cores for which we do detect cm emission, we 
cannot, at this stage, discriminate between free-free emission from young stellar objects 
and spinning dust emission. We emphasise that, although the CARMA observations 
described in this analysis place important constraints on the presence of spinning dust 
in cold, dense environments, the source sample targeted by these observations is not 
statistically representative of the entire population of Galactic cores. 

Key words: radio continuum: ISM - infrared: ISM - ISM: general - ISM: clouds - 
ISM: dust, extinction - ISM: evolution 


1 INTRODUCTION 

In recent years, there has been mounting evidence for the 
existence of a new emission component of the Galactic inter¬ 
stellar medium (ISM). This new component was first iden¬ 
tified in the diffuse ISM, and due to its unexpected discov¬ 
ery, was simply referred to as anomalous microwave emis¬ 


sion (Leitch et al. 1997). We now believe that this anomalous 


microwave emission is due to electric dipole emission from 


* ESA Research Fellow 
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small interstellar dust grains, commonly known as spinning 
dust emission (Dr aine fc Lazarian||1998 ). 

Spinning dust emission occurs in the wavelength 
range ~3 - 0.3 cm and produces a very distinctive peaked 
spectrum, peaking at wavelengths of ~1 cm, with the exact 
wavelength depending on the physical environmental condi¬ 
tions. Although this microwave emission from small spinning 
dust grains has been found in a wide variety of environ¬ 
ments, the best examples of spinning dust detections occur 


in dense environments: the Perseus molecular cloud (Watson 


et al.||2005l |Tibbs et al.||2010[ [2 013 i Pl anck Collaboration 


XX |2011| |Genova-Santos et al.||2015|; and the dark cloud 
LDN1622 ( Finkbeiner|2004 Casassus et al.|2006 Harper et 
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Table 1 . Nummary Q f the 15 cold clumps targeted in this analysis. 


CARMA rms 


Target 

R.A. 

(J2000) 

Decl. 

(J2000) 

Distance 

(kpc) 

Tge. s 

(K) 

^Short 

(mjy beam -1 ) 

^Long 

(mJy beam -1 ) 

ECC181 G102.19+15.24 

20:41:10.74 

+67:21:44.3 

0.33° 

9.6 

0.58 

0.55 

ECC189 G103.71+14.88 

20:53:30.29 

+68:19:32.9 

0.29“ 

9.7 

0.37 

0.37 

ECC190 G103.77+13.90 

21:02:09.19 

+67:45:51.8 

0.29“ 

11.1 

0.45 

0.48 

ECC191 G103.90+13.97 

21:02:23.24 

+67:54:43.3 

0.29“ 

11.1 

0.90 

0.86 

ECC223 G113.42+16.97 

21:59:59.03 

+76:34:08.7 

0.99** 

8.9 

0.49 

0.47 

ECC224 G113.62+15.01 

22:21:37.34 

+75:06:33.5 

0.86 6 

8.0 

0.75 

0.78 

ECC225 G113.75+14.90 

22:24:16.23 

+75:05:01.8 

0.88k 

8.6 

0.67 

0.61 

ECC229 G114.67+14.47 

22:39:35.57 

+75:11:34.0 

0.77k 

10.3 

0.45 

0.38 

ECC276 G127.88+02.66 

01:38:39.14 

+65:05:06.5 

1.16k 

12.6 

0.33 

0.37 

ECC332 G149.41+03.37 

04:17:09.10 

+55:17:39.4 

0.18 6 

8.7 

0.52 

0.59 

ECC334 G149.58+03.45 

04:18:23.96 

+55:13:30.6 

0.20k 

8.7 

0.52 

0.59 

ECC335 G149.65+03.54 

04:19:11.28 

+55:14:44.4 

0.17k 

8.1 

0.57 

0.65 

ECC340 G151.45+03.95 

04:29:56.29 

+54:14:51.7 

0.19“ 

10.1 

0.36 

0.39 

ECC345 G154.07+05.09 

04:47:23.41 

+53:03:31.4 

0.34k 

7.4 

0.59 

0.63 

ECC346 G154.07+05.21 

04:47:57.83 

+53:07:51.2 

0.23k 

10.0 

0.67 

0.64 


“Distances are based on association with known Lynds dark nebulae ( Hilton Sz Lahnlla|l995|>. ^Dista nces are kinematic distances based 
on 13 CO observations ( Wn et al.|2012| . Tg as values taken from |Wu et al. (2012). 


ah||2015 ). In this analysis we attempt, for the first time, to 


search for spinning dust emission in even more dense envi¬ 
ronments: Galactic cores. Galactic cores represent one of the 
earliest phases of star formation, created by the gravitational 
collapse of an over density within a molecular cloud (see 
Bergin & Tafalla 2007, for a complete review). Although no 


such observations of spinning dust emission from Galactic 
cores have been performed, Ysard et al. (2011) conducted 
radiative transfer modelling of such dense cores and pre¬ 
dicted that spinning dust emission should be detectable at 
microwave frequencies if there are small dust grains present. 
Additionally, Ysard et al. (2011) showed that as the den¬ 
sity of the cores increases, an anti-correlation between the 
microwave emission and the mid-infrared (IR) emission oc¬ 
curs. The mid-IR emission is concentrated around the edge 
of the dense cores, while the microwave emission is contained 
within the centre of the cores. Ysard et al. (2011) proposed 
that this anti-correlation could be explained in two ways. Ei¬ 
ther the small grains that produce both the mid-IR emission 
and the microwave emission are present throughout the core, 
but due to the high densities, the stellar photons cannot pen¬ 
etrate deep enough into the core to cause the small grains 
to emit at mid-IR wavelengths, or that there is a deficit 
of small grains in the centre. Therefore, the motivation for 
this current work was to search for spinning dust emission 
in dense cores, and determine which of these two hypothe¬ 
ses is correct. To this end we observed a sample of Galactic 
cores with the CARMA interferometer at a wavelength of 
1 cm, which based on the work by Ysard et al. (2011), is the 
expected peak wavelength of the spinning dust emission in 
these dense environments. 


This paper is organized as follows. In Section[2]we define 
our sample, and in Section [ 3 ] we describe the data that we 
use to perform this analysis. In Section[4]we characterize the 
physical properties of our cores, and in Section[5]we compare 
our observations with the expected level of spinning dust 
emission. Finally, we present our conclusions in Section [6] 


2 SOURCE SELECTION 


For the dust temperatures characteristic of Galactic cold 
cores (~8 - 14 K), the peak of the thermal dust emis¬ 
sion occurs around ~200 /im, implying that the best wave¬ 
lengths to detect these objects is in the far-IR and sub-mm 
bands. Planck , with its all-sky coverage at sub-mm wave¬ 
lengths, is an ideal instrument with which to detect cold 
cores. In fact, as part of the Planck Early Release Com¬ 
pact Source Catalogue ( jPlanck Collaboration VII |2011 ), the 
Early Compact Core (ECC) catalogue was published. The 
ECC catalogue contains 915 cold clumpQ identified by the 
Planck team to be the coldest (T < 14 K) and most re¬ 
liable (SNR > 15) sources in the Cold Core Catalogue of 
Planck Objects (C3PO). 

We selected a sample of 15 sources from the ECC cat¬ 
alogue. The sources were selected based on the following 
criteria: (a) the source must have been observed with the 
Herschel Space Observatory PACS and SPIRE photome¬ 
ters (Pilbratt et al. 2010); (b) the source must be visible 
from the site of the CARMA interferometer; (c) the source 
must have an angular size less than 12 arcmin to ensure that 
we do not resolve out any flux; (d) the source must have a 
known distance estimate; and (e) the sources must cover a 
range of densities, nu- The 15 sources included in our sample 
are listed in Tabled] 

Also listed in Tableware the distance estimates for each 
source. [Wu et al.| (2012) observed the J — 1 — 0 transition 
of 12 CO, ii5 CO, and C^O for 674 of the 915 ECC sources 
using the 13.7 m telescope of the Purple Mountain Observa¬ 
tory. Using the 13 CO Ulsr measurements with the Clemens 


(1985) rotation curve, and adopt i ng va lues of Ro = 8.5 kpc 


and Go = 220 km s , Wu et al. (2012) computed the kine¬ 


matic distance to each source. For sources within the Solar 


1 Since the angular resolution of Planck is ~5 arcmin in the far- 
IR/sub-mm bands, it can not detect individual cores, and there¬ 
fore we refer to the Planck sources as clumps. 
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circle, which suffer from the kinematic distance ambiguity, 
they selected the near distance. Since Wu et al. (2012) was 
unable to assign a distance estimate to all of their ECC 
sources, there were five sources in our sample (ECC 181, 
ECC189, ECC190, ECC191, and ECC340) for which we did 
not have a kinematic distance estimate, and for these sources 
we estimated the distance based on associations with known 
Lynds dark nebulae (Hilton & L ahulla||1995 ). 

In addition to a distance estimate, Wu et al. (2012) 
also derived the temperature of the gas, T gas , for all of their 
ECC sources. The T gas values for our 15 targets are also 
listed in Table [l] and these values are adopted throughout 
this analysis. 


3 OBSERVATIONS 
3.1 CARMA Data 

CARMA is an interferometer consisting of a total of 23 an¬ 
tennas located at a high-altitude site in eastern California. 
CARMA operates at three bands: 1 cm (27 - 35 GHz), 
3 mm (85 - 116 GHz), and 1 mm (215 - 270 GHz), and 
consists of six 10.4 m, nine 6.1 m, and eight 3.5 m an¬ 
tennas that are combined to form the full CARMA array. 
The data presented in this analysis were observed in the 
1 cm band with fifteen 500 MHz channels between 27.5 and 
35 GHz using the eight 3.5 m antennas. Six of the eight an¬ 
tennas formed a compact configuration with baseline lengths 
of ~0.3 - 2.0 kA, while the remaining two antennas pro¬ 
vided much longer baselines of ~2.0 - 8.0 kA. We refer to 
data from the baselines between the six compact antennas 
as short-baseline data, and data from the baselines between 
the six compact antennas and the two outlying antennas as 
long-baseline data. In Figure [l] we plot the Fourier space u,v 
coverage for our observations of one of our targets, ECC181. 
The u,v coverage of the rest of our fifteen targets is very sim¬ 
ilar. In the left panel of Figure [l] we plot the complete u,v 
coverage for all baselines, which clearly illustrates how we 
separate the data into short- and long-baseline data, while 
in the right panel we plot a zoomed in view of the short 
baselines. 

Our CARMA observations were performed during three 
separate observing runs in May 2012, August 2012, and Jan¬ 
uary 2013. The data were reduced using a suite of MATLAB 
routines, which convert the data to physical units, correct 
for instrumental phase and amplitude variations, flag data 
that do not meet designed criteria at each stage of reduction, 
and perform the data calibration (see |Muchovej et al.|2007 
for further details on the reduction pipeline). The absolute 
calibration is derived from observations of Mars, and based 
on long-term monitoring of flux density calibrators, we esti¬ 
mate that the calibration is accurate to 5%. The reduction 
pipeline outputs the calibrated visibility data, which were 
then imported into AIPS, where the task IMAGR was used 
to perform the Fourier transform and deconvolution. This 
task uses a CLEAN-based algorithm to perform the decon¬ 
volution, and we implemented a Briggs robust parameter 
of 0, providing a compromise between uniform and natural 
weighting. We imaged the short-baseline data and the long- 
baseline data separately, using a 120 arcsec and 20 arcsec 
restoring beam (FWHM), respectively. The final CARMA 


maps are displayed in Figure [AT] and the rms noise in both 
the short-baseline (crshort) and long-baseline (crLong) images 
for each source is tabulated in Tabled] 

The long-baseline data are ideal for identifying com¬ 
pact radio sources (< 1 arcmin), while the short-baseline 
data are sensitive to the more diffuse emission on angular 
scales of ~2 - 12 arcmin. Looking at the maps in Figure [Al] 
we detected four compact sources that are visible in both 
the short- and long-baseline data (one source each in the 
ECC189, ECC191, ECC340, and ECC346 maps). Using the 
AIPS task JMFIT we fitted 2-dimensional Gaussians to each 
source to derive positions and flux densities. We searched 
the literature to identify these sources, and the NVSS coun¬ 
terpart for each source is listed in Table [2] Also listed in 
Tableware the measured CARMA 1 cm flux densities. The 
CARMA flux density uncertainties include a 5% calibration 
uncertainty combined in quadrature with the uncertainty 
from the fit to the source. Throughout the rest of this anal¬ 
ysis we ignore these four point sources. 

In addition to the point sources, there is not much ex¬ 
tended emission present in the short-baseline maps displayed 
in Figure [AT] In fact, we find that only 1 of our 15 clumps, 
ECC229, has a significant detection of emission, which con¬ 
sists of a more compact component to the west and a more 
extended structure, likely consisting of two separate compo¬ 
nents, to the east. 


3.2 Herschel Data 


In order to characterize the physical conditions of the ob¬ 
served CARMA clumps, we used Herschel PACS and SPIRE 
data. Eight of our targets are covered by the Guaranteed 
Time Key Programme “Probing the origin of the stellar ini¬ 
tial mass function: A wide-held Herschel photometric survey 
of nearby star-forming cloud complexes” (KPGT_pandre_l; 
Andre et al. 2010) while the remaining seven are covered 


by the Open Time Key Programme “Galactic Cold Cores: 
A Herschel survey of the source populations revealed by 
Planck” (KPOT_mjuvela_l; Juvela et al.| 2010). 

The Herschel data were downloaded from the Herschel 
Data Archive and re-processed using the Herschel Interac¬ 
tive Processing Environment (HIPE) version 10.0.2747. The 
PACS data were processed using the standard pipeline for 
extended sources. Starting with the level 1 PACS data, this 
processing performed a correction for the global correlated 
signal drift and produced the time ordered data, which were 
used to create a map using the Madmap map-maker. Finally, 
a point source artifact correction was applied, resulting in 
maps at 70, 100, and 160 /xm with an angular resolution 
of-6, 7, and 11 arcsec, respectively. The absolute calibration 
uncertainty of the PACS data was conservatively estimated 
to be 10 %. The SPIRE data were processed from level 1 
to level 2.5 using the SPIA reduction routines, with the De¬ 
striper map-maker used to produce maps at 250, 350, and 
500 fi m with an angular resolution of ~18, 25, and 36 arc¬ 
sec, respectively. The absolute calibration uncertainty of the 
SPIRE data was conservatively estimated to be 7 %. 

Although PACS has three channels (70, 100, and 
160 (i m), only two (either 70 or 100 /im along with 160 /xm) 
can be observed simultaneously for each observation. There¬ 
fore, to ensure that observations in all three PACS channels 
are obtained, a few of our sources had both PACS/SPIRE 
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u -z -o ^.U 1.0 l.u U.O U.U -U.o -1,1 

u (Kilo Wavelength) u (Kilo Wavelength) 


Figure 1 . Complete u-v coverage (left) and the u-v coverage for only the short-baselines (right) for our observations of ECC181. The 
colour scale represents the 15 channels with red corresponding to the highest frequency channel and blue corresponding to the lowest 
frequency channel. The u,v coverage for the rest of our targets is very similar. 


Table 2. Flux densities of the four point sources identified in the CARMA maps. 


Source 


Short Baselines 


Long Baselines 

R.A. 

(J2000) 

Decl. 

(J2000) 

Si cm 

(mjy) 

R.A. 

(J2000) 

Decl. 

(J2000) 

Si cm 

(mjy) 

NVSS 205303+682200 

20:53:04.00 

+68:22:06.4 

2.85 ± 0.64 

20:53:04.44 

+68:21:58.0 

3.45 ± 0.65 

NVSS 210243+675819 

21:02:43.91 

+67:58:16.7 

101.77 ± 5.30 

21:02:43.70 

+67:58:17.1 

102.52 ± 5.27 

NVSS 043015+541429 

04:30:15.40 

+54:14:31.6 

4.41 ± 0.62 

04:30:14.84 

+54:14:28.1 

4.32 ± 0.69 

NVSS 044819+530830 

04:48:18.88 

+53:08:22.9 

10.26 ± 1.40 

04:48:19.47 

+53:08:30.4 

9.98 ± 1.07 


parallel observations and PACS photometry observations. 
However, this not only results in observations in all three 
PACS channels, it also results in a duplication of observa¬ 
tions in the 160 /im channel. Not all of our sources had du¬ 
plicate 160 /im observations, but for the ones that did, we 
used the mean of the two 160 /im maps. All fifteen sources 
did have at least one set of PACS and SPIRE observations, 
meaning that we have Herschel observations at 160, 250, 
350, and 500 /im for all of our sources which, as we will 
discuss in Section [471] is adequate for our analysis. 


4 CLUMP PROPERTIES 

To understand why only 1 of our 15 clumps had detectable 
1 cm emission, it is crucial to explore the physical conditions 
of our sample. Properties such as the density and radiation 
field are important parameters in spinning up the grains 
to produce spinning dust emission. Additionally, we need 
to determine the evolutionary stage of the clumps, specifi¬ 
cally, whether these sources are gravitationally bound and, 
if so, whether they are pre-stellar or proto-stellar. If indeed 
some of the clumps are proto-stellar, we also need to iden¬ 
tify the location of young stellar objects (YSOs) as they 
can generate emission at cm wavelengths due to a variety 
of mechanisms such as stellar winds and/or shock-induced 
ionization, which can mimic the spinning dust emission spec¬ 


tral signal. Additionally, knowledge about the evolutionary 
stage of the clumps and, in particular, understanding if the 
observed sources are both simultaneously forming stars and 
harbouring spinning dust emission can help to understand 
the potential role of spinning dust emission in the star for¬ 
mation process. 

Therefore, in Section |4Tl we compute the hydrogen col¬ 
umn density and dust temperature within our clumps, and 
in Sect ion [472] we estimate the hydrogen number density and 
the radiation field. Then in Section 14731 we estimate the mass 
and determine if the clump is stable or likely to undergo col¬ 
lapse, and in Section [4.4| we identify and classify candidate 
YSOs in the vicinity of our clumps. 


4.1 Hydrogen Column Density and Dust 
Temperature 

All of the Herschel maps were convolved to the angular res¬ 
olution of the 500 /im map using the convolution kernels 


produced by Aniano et al. (2011) to allow for accurate com¬ 


parisons between each of the Herschel bands. To ensure we 
are probing the dense cores themselves and are not contam¬ 
inated by foreground or background emission, we performed 
a background subtraction on the Herschel maps. For each 
of the Herschel maps, we computed the median value of the 
flux in a reference position that is devoid of emission, and 
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subtracted this from the map. The emission in these back¬ 
ground subtracted Herschel maps was then modelled using 


I u — fimnNKKjsBviTd), 


(i) 


where I v is the intensity at frequency v, /x is the molecular 
weight per hydrogen atom, which we assume to be 1.4, ran is 
the mass of a H atom, Ah is the hydrogen column density, 
k v is the dust opacity, and B u (Td) is the Planck function 
for dust temperature Td. The dust opacity is the subject of 
much debate, with large variations observed between mod¬ 
els, depending on the physical properties of the dust grains 


such as size, composition, and structure (e.g., Ossenkopf & 
Henning 1994). In this analysis we use the dust opacity pa¬ 
rameterization defined by Beckwith et al. (1990) 


Av i/ — 0 


.if_r_V 

V1000 GHz/ 


2 -1 
cm g , 


( 2 ) 


which is applicable to these cold, dense environments (e.g., 
Ward-Thompson et al.||2010| [Planck Collaboration XXII| 


2011||Planck Collaboration XXIII|2011| jJuvela et al.|2012| ). 

This normalization assumes a standard gas to dust mass 
ratio of 100. 

We used the model defined by Equations 0 and 0 to 
fit the background subtracted Herschel maps at 160, 250, 
350, and 500 /xm. We excluded the 70 and 100 /xm Herschel 
maps from the fit because our clumps are cold (< 14 K; this 
was one of the selection criteria used to produce the Planck 
ECC catalogue), and therefore they do not emit strongly 
at wavelengths ^ 100 /xm. Furthermore, excluding the 70 
and 100 /xm maps ensures that we avoid possible contami¬ 
nation from stochastically heated small dust grains. We in¬ 
spected the 70 and 100 /xm maps for all of the clumps and 
confirmed that there is very little emission present at these 
wavelengths, justifying our decision to exclude them from 
the fit. 

We fitted the data for Ah and Td , fixing the dust opacity 
index (3 = 2 (e.g., Planck Collaboration XXII||2011 Planck 


Collaboration XXilI||2011| |Juvela et al.||2012|). This fit 


performed on a pixel-by-pixel basis using MPFIT ( Mark- 
war dt| 2009 ) resulting in maps of Ah and Td, along with their 
associated uncertainty maps, which were computed from the 
covariance matrix. The Ah and Td maps for all 15 clumps 
are displayed in Figure [Bl] For our 15 targets, we find val¬ 


ues of Ah ranging from ~6xl0 to ~7xl0 H cm and 
dust temperatures between Td ~10 to ~15 K. These val¬ 
ues are well within the the range found for the entire C3PO 


catalogue by Planck Collaboration XXIII (2011), and are 


consistent with those estimated by Juvela et al. (2012) us¬ 


ing Herschel SPIRE data for their sample of cores. 


4.2 Hydrogen Number Density and Radiation 
Field 

Herschel , with its sub-arcmin angular resolution, is able to 
resolve the sources observed by Planck into sub-structures. 
In the case of the ECC clumps, we can therefore identify 
individual cores within each clump. For this purpose, we 
used the clumpfind package ( Williams et al.| [l994), which 
identifies local peaks in the Ah maps and follows them down 
to lower levels, resulting in a decomposition of the map into 


multiple cores, clumpfind identified 34 cores, the positions 
of which are listed in Table [ 3 ] and labelled in Figure |B1| 
In addition to the position, CLUMPFIND also computes the 
angular size of each identified core. Since this angular size 
is the effective circular radius, R (= (A p i X eis/7r) 1//2 , where 
-Xpixeis is the total number of pixels), we assumed that the 
cores are spherical, and for each core we computed the mean 
column density, Ah, and the mean dust temperature, Td- 
The values of Ah and Td for each core are listed in Table [3] 
Using the distance to each clump (see Table [l]) we con¬ 
verted the angular size of each core computed by CLUMP 
FIND to a physical, linear size, L. Since we have performed a 
background/foreground subtraction, by combining Ah and 
the linear size we can estimate the mean density, txh- Based 
on the conservation of mass, the total mass in the Ah map 
within a given radius, T, is equivalent to the mass within a 
sphere of the same radius, with a mean density, txh, he., 


7tT 2 /X?7Xh^Vh 


4 

3 


irR 3 /xttxh^h- 


Since R = L/2, this results in 


(3) 


3 Nh 

Uh= 2UT- 


(4) 


Using Equation 0 we estimated txh for each core. 
The uncertainty on nu was calculated assuming a 10% 
uncertainty on the distance. The derived txh values 
range from ~5xl0 3 H cm -3 in core 1 in ECC223 
to ~125xl0 3 H cm -3 in core 2 in ECC229 (see Table [3]). 
Although these values are mean quantities, they are com¬ 
parable to the densities estimated in the entire C3PO cata¬ 
logue (Planck Collaboration XXlil 2011). 

In addition to estimating txh, we also estimated the 
mean radiation field, Go, in each core by converting Td into 
Go using 


Go = 


17.5 K ) 


(5) 


where we fixed = 2 as before. The computed values of 
Go are also listed in Table [3] and we find values rang¬ 
ing from ~0.04 to 0.30, where a value of 1 corresponds to 


the Mathis et al. (1983) solar neighbourhood radiation field. 


Such low values of Go indicate that the interior of these 
dense environments are shielded from the surrounding in¬ 
terstellar radiation field. 


4.3 Mass 

We estimated the total (dust + gas) mass for each of our 
cores using Equation 0- The computed mass estimates are 
listed in Tableland range between ~0.4 - 115 M 0 . 

To determine if our cores are gravitationally stable or 
in the process of collapsing, we also computed the Bonnor- 
Ebert mass, Mbe, the mass of an isothermal sphere in hy¬ 
drostatic and pressure equilibrium. Cores with M > Mbe 
are unstable, and will therefore collapse, while cores with 
M < Mbe are stable. We estimated the Bonnor-Ebert mass 
as 
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Mbe ~ 1 


.82 ( 


n H 

10 4 H cm- 


10 K 


Mr; 


O’ 


(6) 


from Lada et al. (2008). Both the values of M and Mbe for 
all of the cores are listed in Table |3] We find that 29 of our 
34 cores are unstable, and therefore undergoing collapse. 


4.4 Young Stellar Objects 

Since the unstable cores will collapse and eventually form 
stars, we also investigated the ongoing star formation by 
searching for YSOs in the vicinity of our cores. Various IR 
colour-cuts have been used to identify YSOs, e.g., based on 
Spitzer (e.g., Allen et ah1|2004 ) and WISE (e.g., Koenig et 
al. 2012) data. Since not all of our clumps have been ob¬ 


served with Spitzer , we take advantage of the all-sky cover¬ 
age of WISE, which mapped the entire sky simultaneously 


at four wavelengths: 3.4, 4.6, 12, and 22 /xm (Wright et al. 


2010). We apply the colour-colour selection criteria devel¬ 


oped by Koenig et al. (2012) to the AllWISE all-sky source 
catalog to identify YSO candidates, excluding sources with 
a signal-to-noise ratio < 5 and contaminated sources iden¬ 
tified by the catalog contamination flags “D” (diffraction 
spikes), “P” (persistent latent artefacts), “H” (halos from 
bright sources), and “O” (optical ghosts from bright sources). 
This selection criteria allows us to reject extragalactic con¬ 
taminants and distinguish between class I and class II YSO 
candidates (for details of the specific selection criteria, see 
Appendices A.l, A.2, and A.3 in Koenig et al.|2012 ). 

This analysis enabled us to identify YSO candidates as¬ 
sociated with each clump. We found that only 5 of our 15 
clumps have at least one YSO candidate within our 20 ar- 
cmin x 20 arcmin maps displayed in Figure |B1[ and that 
only 4 of these (ECC190, ECC191, ECC229, and ECC276) 
have a YSO candidate within the CARMA primary beam. 
The location of the YSO candidates, and their classifica¬ 


tion (either class I or class II), are marked on Figure B1 It 
is widely accepted that no colour-colour identification crite¬ 
ria are 100% accurate. Using the contamination rate of false 


identifications from Koenig et al. (2012), we estimate that 


there are ~0.4 - 0.9 false identifications in each of our maps. 

Looking at the location of these YSO candidates rela¬ 
tive to our cores, we determined that 6 cores are possibly 
associated with a YSO (see Table [3]). Core 1 in ECC191 is 
coincident with a class I YSO candidate, while cores 2, 3, 
and 4 in ECC229 are coincident with two class I YSO candi¬ 
dates, a class II YSO candidate, and a class I YSO candidate, 
respectively. In addition, there are 2 cores which have YSO 
candidates nearby: core 1 in ECC229 has two nearby class 
ii YSO candidates, and core 2 in ECC276 has a class I YSO 
candidate nearby. 

YSOs have previously been detected at cm wave¬ 
lengths (e.g., AMI Consortiu m et al.| 2011), but of the 6 
cores that are possibly associated with a YSO candidate, 
only three have detectable cm emission (ECC229 cores 1, 2, 
and 3). However, given the variety of mechanisms by which 
YSOs can produce cm emission, perhaps it is not surprising 
that we do not detect cm emission from all of our candi¬ 
dates. This also makes it more difficult to interpret the cm 
emission that we detect in ECC229, as we will discuss in 
Section 15.11 


5 DISCUSSION 


5.1 ECC229 

Out of all of our 15 clumps, ECC229 is the only clump for 
which we detect any extended cm emission, making it par¬ 
ticularly interesting. In our CARMA map of ECC229 we 
detected three components of cm emission (associated with 
cores 1, 2, and 3), although we identified four cores. Based 
on the values listed in Table [3] the cores in ECC229 have 


B1 


the highest values of N h- Moreover, as shown in Figure 
the cm emission appears to be spatially correlated with Affi 
and originating from the densest regions. 

Core 2 in ECC229 has the highest density (% = (123.7 
zb 17.3) x 10 3 H cm -3 ), but also the warmest dust tem¬ 
perature (?d = 14.3 =b 0.3 K) and highest radiation field 
value (Go = 0.298 zb 0.041), potentially suggesting that the 
star formation process has already begun. Indeed two YSO 
candidates, likely heating the dust in the core, are found to 
be associated with core 2. In fact all 4 of the cores in ECC229 
appear to be associated with a YSO candidate, and the fact 
that we only detect cm emission from 3 of the cores might 
suggest that the cm emission is not linked to the presence of 
the YSOs. However, given that YSOs can produce cm emis¬ 
sion via a variety of different mechanisms (accretion, stellar 
winds, shock ionization), it is far from clear. For example, 
the spectral index at cm wavelengths (defined as a , where 
St, oc predicted from free-free emission from an ionized 
stellar wind can cover the range from —0.1 ^ a ^ +2.0, with 
typical values of a ~0.6 in the case of a spherical, isother¬ 


mal stellar wind (Reynolds 1986). A similar range of val¬ 
ues for a are produced for free-free emission from accretion 
and shock ionization. The spectral index from spinning dust 
emission at these wavelengths is dependent on the frequency 
of the peak of the emission, and so a can be both positive or 
negative. This makes it difficult to distinguish between the 
cm emission arising from YSOs as opposed to spinning dust 
emission, especially in the absence of multi-wavelength cm 
observations, which can closely trace the microwave emis¬ 
sion, combined with IR spectroscopic data, which can pro¬ 
vide details on the YSO geometry (e.g., disk inclination, 
etc.). Therefore, lacking this ancillary information, to de¬ 
termine the presence of spinning dust emission (or a lack 
thereof), we compared our CARMA observations to the pre¬ 
dicted level of spinning dust emission based on the physical 
conditions in each core. 


5.2 Spinning Dust Emission 

With the physical properties of the clumps characterized, as 
described in Section [4] we can investigate the significance 
of our CARMA observations, namely: does the fact that we 
did not detect cm emission in 14 of our 15 clumps signify 
that there is no spinning dust emission present? In order 
to address this question, we computed the expected level of 
spinning dust emission, using the constraints provided by 
the analysis of the Herschel data. 

We characterized our cores by using the hu and Go val¬ 
ues derived in Section |+2| the gas temperatures, T gas , listed 
in Table [I] and the standard values for dense environments 
for the hydrogen ionization fraction (xu = 0), the carbon 
ionization fraction (xc = 10 -6 ), and the molecular hydrogen 
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Figure 2. Spectrum for the four cores in ECC229, consisting of both a spinning dust ( dot-dashed line) and thermal dust (dotted line) 
component. The thermal dust emission is modelled using a modified black body with /3 fixed at 2, and is fitted to the integrated flux 
density at 160, 250, 350, and 500 jum ( diamonds) in each core, while the spinning dust curve, and associated uncertainty (shaded region), 
is the predicted level of spinning dust emission based on the modelling analysis. Also plotted is the measured flux density at 1 cm (for 
core 4 in which we do not detect any cm emission, we have plotted the 5cr upper limit), illustrating how the predicted spinning dust 
emission is larger than what is observed with CARMA. 


fraction (y — 0.999) as defined by Draine & Lazarian (1998). 
In addition to these parameters, and since these are dense 
environments, we used the Weingartner & Draine (2001) 
grain size distribution for R\ = 5.5 with be = 3 x 10 5 (cor¬ 
responding to the parameters in [Weingartner &; Draine|2001[ 
Table 1, Line 16). With these parameters as inputs, we used 
the spinning dust model, SPDUST (Silsbee et al. 2011) to pre¬ 
dict the expected level of spinning dust emission at 1 cm. 


SPDUST estimates the spinning dust emissivity in units 
of Jy sr -1 H -1 cm 2 which we converted to a flux by using 
our computed value of Nn. To estimate the uncertainty on 
the predicted spinning dust emission, we repeated this mod¬ 
elling analysis 1000 times, randomising the input nn and 
Go values within their uncertainty, resulting in a range of 
predicted values for the spinning dust emission at 1 cm. We 
then compared the predicted spinning dust emission with 
our CARMA maps. For the cores for which we do not de¬ 
tect any cm emission we computed a conservative 5cr upper 
limit by scaling the rms noise in the CARMA maps based on 
the ratio of the size of the beam to the size of the core, while 
for ECC229, we fitted the extended emission with three 2- 
dimensional Gaussian components using the AIPS task JM- 
fit. The last two columns in Table [3] list our measured 1 cm 
flux densities (5'i b c S m Ved ) along with the predicted levels of 
the spinning dust emission (Sf r ^ cted ). For all of our cores 


we find that the predicted level of spinning dust emission is 
much larger than what we observe in our CARMA maps. 
We illustrate this result for the four cores in ECC229 in Fig¬ 
ure^ where we plot the predicted spinning dust emission, 
with its associated uncertainty, the thermal dust emission, 
modelled using a modified black body with /3 fixed at 2, and 
our CARMA 1 cm flux density measurements. 


We do note, however, that this result depends strongly 
on the spinning dust modelling analysis. Firstly, since we are 
only comparing our observed cm emission with the predicted 
spinning dust emission, we are actually being conservative 
as we are ignoring any possible contribution to the 1 cm 
emission from thermal dust emission or free-free emission, 
which would increase the expected level of the 1 cm emis¬ 
sion. Secondly, it is possible that the discrepancy between 
the observed and the predicted level of spinning dust emis¬ 
sion could be due to an incorrect assumption in our spinning 
dust modelling. For example, there is much uncertainty on 
the value of the electric dipole moment of dust grains. The 
value used in this analysis is consistent with the standard 
value estimated by Draine & Lazarian (1998) - an aver¬ 
age rms dipole moment per atom of 0.38 D. If this value 
is overestimated, it would result in an overestimate of the 
level of the spinning dust emission. Likewise, we might have 
overestimated as a consequence of the uncertainty in the 
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distance estimate or due to our assumption that the core ge¬ 
ometry is spherical. Similarly, our estimate of Th, and hence 
Go , may be overestimated due to the mixing of different line- 
of-sight temperature variations (e.g. Nielbock et al. 2012 
Roy et al. 2014). Again, such overestimates would result 


us overestimating the spinning dust emission. Finally, our 
adopted grain size distribution might not be appropriate for 
Galactic cores. Although we used the grain size distribution 
for dense environments with Ry = 5.5 (as opposed to values 
of Ry = 3.1 or Ry = 4.0), this could still be overestimat¬ 
ing the abundance of the small grains responsible for the 


spinning dust emission. In a forthcoming paper (Tibbs et al. 


2015) we will investigate these scenarios in more detail and, 


in particular, use these observations to place a constraint 
on the abundance of the small dust grains in our sample of 
cores. 

Finally, we stress that the sample of Galactic cores 
which was the subject of this study does not span the full 
parameter space (in terms of density, temperature, mass, 
etc.) relative to the Galactic cold core population. For this 
reason, the conclusions reached by our analysis cannot eas¬ 
ily be generalized, and more cm observations like the ones 
presented are needed to allow confirmation of our findings. 
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6 CONCLUSIONS 

We have attempted, for the first time, to search for spinning 
dust emission in a sample of Galactic cold clumps. To do 
this we observed a sample of 15 cold clumps with CARMA 
at 1 cm, and found that only 1 of our 15 clumps exhibited 
significant extended emission. To determine if the lack of de¬ 
tection of 1 cm emission could rule out a spinning dust detec¬ 
tion, we investigated the physical properties of our sample 
of clumps using ancillary Herschel and WISE data. Using 
Herschel photometric data at 160, 250, 350, and 500 (i m of 
our 15 cold clumps we produced maps of Nu and Xd, from 
which we identified 34 cores, characterized by the densest 
and coldest regions within each of the 15 clumps. For each 
core we estimated Ah, Td, mp and Go. Making use of the 
all-sky coverage of WISE, we used the AllWISE source cat¬ 
alog to identify candidate YSOs and found that 6 of our 34 
cores were associated with a YSO candidate. 

With the physical environments of each core con¬ 
strained, we used SPDUST to model the spinning dust emis¬ 
sion, which we compared to our CARMA observations, and 
we found that the observed cm emission in all of our cores 
was below the predicted level. This implies that we do not 
detect spinning dust emission from 14 of our 15 clumps, and 
in the 1 clump in which we do detect cm emission, it could 
be due to either spinning dust emission, but at a level much 
lower than predicted based on the modelling, or it could be 
due to free-free emission from YSOs. 

This analysis is the first attempt to detect spinning dust 
emission in Galactic cores, however, we emphasize that our 
sample is not statistically representative of the entire Galac¬ 
tic cold core population, and therefore we are cautious to 
extend this result to all of the cold cores in the Galaxy. 
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Table 3. Derived properties of the cores analysed in this work, including the core size, mean column density (A/’h), mean density (nn), mean dust temperature (T^), mean radiation 
field (Go), mass (M), Bonnor-Ebert mass ( Mbe ), associated YSOs (YSO), observed 1 cm flux density (£i b c S m rved )> an d the predicted 1 cm spinning dust emission (5^ r ^ lcted ). 


Clump 

Core 

R.A. 

(J2000) 

Decl. 

(J2000) 

Size 

(pc) 

Ah 

(10 2o H cm“ 2 ) 

nn 

(10 3 H cm -3 ) 

T d 

(K) 


Go 


M 

(M©) 

Mbe 

(M 0 ) 

YSO 

cobserved 
°1 cm 

(mJy) 

^predicted 
^1 cm 

(mJy) 


1 

20:41:13.2 

67:20:34.2 

0.23 

121.1 zb 10.5 

25.6 dz 3.4 

11.6 dz 0.2 

0.086 

zb 

0.008 

5.6 zb 0.9 

1.1 

_ 

< 2.87 

9.09 - 

- 12.16 

ECC181 

2 

20:40:56.2 

67:22:54.4 

0.18 

115.2 db 10.0 

32.0 dz 4.2 

11.7 db 0.2 

0.090 

zb 

0.009 

3.1 zb 0.5 

1.0 

- 

< 1.66 

5.10 

- 6.79 


3 

20:40:31.9 

67:20:48.1 

0.13 

107.4 ± 9.3 

41.3 dz 5.5 

11.4 zb 0.2 

0.078 

zb 

0.007 

1.5 zb 0.2 

0.8 

- 

< 0.86 

2.53 

- 3.34 

ECC189 

1 

20:53:35.4 

68:19:19.9 

0.17 

65.0 ± 5.7 

18.7 zb 2.5 

12.3 =b 0.2 

0.120 

zb 

0.012 

1.6 zb 0.3 

1.3 

- 

< 1.30 

3.47 

- 4.84 

ECC190 

1 

21:01:54.4 

67:43:45.7 

0.20 

128.3 ± 11.1 

31.6 zb 4.2 

12.3 zb 0.2 

0.122 

zb 

0.012 

4.4 zb 0.7 

1.2 

_ 

< 2.15 

10.53 

- 14.04 

2 

21:02:21.5 

67:45:37.7 

0.21 

110.7 ± 9.7 

25.4 dz 3.4 

12.7 =b 0.2 

0.146 

zb 

0.015 

4.4 zb 0.7 

1.3 

- 

< 2.48 

10.25 

- 13.72 

ECC191 

1 

21:02:23.2 

67:54:15.2 

0.21 

274.9 dz 23.8 

63.8 zb 8.4 

11.6 zb 0.2 

0.083 

zb 

0.008 

10.6 zb 1.8 

0.8 

Class I 

< 4.80 

26.62 

- 34.99 

ECC223 

1 

21:59:38.8 

76:33:13.1 

1.06 

115.2 ± 10.2 

5.3 zb 0.7 

12.1 =b 0.2 

0.108 

zb 

0.011 

113.6 dz 18.9 

2.1 

- 

< 5.69 

25.15 

- 48.26 


1 

22:21:27.6 

75:04:26.0 

0.48 

130.7 ± 11.6 

13.4 dz 1.8 

12.2 zb 0.2 

0.112 

zb 

0.011 

25.9 zb 4.3 

1.1 

_ 

< 2.36 

6.29 - 

- 12.85 

ECC224 

2 

22:21:41.2 

75:06:06.0 

0.51 

117.2 ± 10.4 

11.2 ± 1.5 

12.4 zb 0.2 

0.129 

zb 

0.013 

26.8 zb 4.5 

1.2 

- 

< 2.72 

9.59 - 

- 13.79 


3 

22:21:50.5 

75:09:09.5 

0.28 

83.7 ± 7.4 

14.6 ± 1.9 

12.6 db 0.2 

0.138 

zb 

0.014 

5.7 zb 1.0 

1.1 

- 

< 0.81 

1.65 

- 3.13 

ECC225 

1 

22:24:09.1 

75:04:33.1 

0.71 

110.5 zb 9.6 

7.5 zb 1.0 

12.0 db 0.2 

0.105 

zb 

0.010 

49.4 dz 8.2 

1.7 

- 

< 4.49 

14.97 

- 22.46 


1 

22:39:39.5 

75:12:01.5 

0.34 

676.0 ± 57.6 

97.0 zb 12.7 

10.2 db 0.1 

0.040 

zb 

0.003 

68.1 zb 11.3 

0.6 

2 x Class ii 

2.87 zb 0.88 

23.16 

- 30.30 

ECC229 

2 

22:38:47.9 

75:11:27.2 

0.24 

619.5 ± 60.5 

123.7 zb 17.3 

14.3 zb 0.3 

0.298 

zb 

0.041 

32.3 zb 5.6 

0.5 

2 x Class i 

4.12 zb 0.91 

10.83 

- 14.45 

3 

22:39:31.6 

75:11:06.6 

0.35 

642.6 ± 54.8 

89.5 =b 11.8 

10.1 db 0.1 

0.036 

zb 

0.003 

68.8 zb 11.4 

0.6 

Class ii 

2.59 zb 0.96 

23.36 

- 30.59 


4 

22:39:06.5 

75:11:52.5 

0.32 

540.3 ± 47.4 

82.7 =b 11.0 

11.2 =b 0.2 

0.068 

zb 

0.006 

48.0 zb 8.0 

0.7 

Class I 

< 0.78 

16.18 

- 21.28 

ECC276 

1 

01:38:34.7 

65:05:48.7 

0.47 

161.4 db 14.4 

16.7 dz 2.2 

11.1 db 0.2 

0.065 

zb 

0.006 

31.5 zb 5.3 

2.0 

_ 

< 0.56 

4.94 

- 6.70 

2 

01:38:30.3 

65:04:52.7 

0.45 

173.6 ± 15.9 

18.6 dz 2.5 

11.6 zb 0.2 

0.086 

zb 

0.009 

31.4 dz 5.3 

1.9 

Class I 

< 0.52 

4.97 

- 6.74 


1 

04:17:23.9 

55:16:15.7 

0.20 

79.5 ± 7.3 

19.8 zb 2.7 

13.4 zb 0.3 

0.203 

zb 

0.024 

2.7 zb 0.5 

1.1 

_ 

< 6.24 

20.90 

- 37.23 

ECC332 

2 

04:17:04.3 

55:13:55.0 

0.08 

74.5 ± 7.0 

47.2 =b 6.5 

13.4 db 0.3 

0.204 

zb 

0.025 

0.4 zb 0.1 

0.7 

- 

< 0.96 

1.97 

- 2.64 

3 

04:17:02.6 

55:15:47.0 

0.15 

68.2 ± 6.4 

22.1 =b 3.0 

13.6 ± 0.3 

0.218 

zb 

0.027 

1.4 zb 0.2 

1.0 

- 

< 3.69 

10.20 

- 19.20 


4 

04:16:57.5 

55:20:13.2 

0.16 

63.8 zb 5.9 

19.2 =b 2.6 

13.6 zb 0.3 

0.222 

zb 

0.027 

1.5 zb 0.2 

1.1 

- 

< 4.25 

12.86 

- 20.27 

ECC334 

1 

04:18:25.6 

55:12:48.6 

0.29 

89.4 dz 8.2 

14.8 =b 2.0 

13.3 zb 0.3 

0.196 

zb 

0.023 

6.7 zb 1.1 

1.2 

- 

< 11.35 

50.94 

- 72.20 


1 

04:18:51.6 

55:14:44.9 

0.15 

77.7 ± 7.0 

24.4 =b 3.3 

13.3 zb 0.3 

0.191 

zb 

0.022 

1.6 zb 0.3 

0.8 

_ 

< 4.80 

10.43 

- 21.45 

ECC335 

2 

04:19:16.2 

55:14:02.2 

0.16 

91.9 =b 8.2 

28.1 zb 3.8 
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Figure Al. CARMA 1 cm maps for the short (left) and long (right) baseline data for all fifteen cold clumps. The location of the cores 
are numbered and identified by red circles. CARMA contours, starting at 5cr and increasing linearly, are displayed for the short baseline 
data (except for ECC191 where the contours increase in steps of 10cr) along with the CARMA 11 arcmin primary beam (green circle). 
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Figure Al. Continued 
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Figure Al. Continued 
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Figure Al. Continued 
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Figure Al. Continued 
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APPENDIX B: N H AND 7b MAPS 
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Figure Bl. Ah {left) and T ^ {right) maps for all fifteen cold clumps. The location of the cores are numbered and identified by red 
circles, while class I and class II YSO candidates are labelled as blue and green stars, respectively. CARMA short-baseline data contours, 
starting at 5cr and increasing linearly, are overplotted (except for ECC191 where the contours increase in steps of 10cr) along with the 
CARMA 11 arcmin primary beam {green circle). 
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Figure Bl. Continued 


© RAS, MNRAS 000, ?? [TT] 














































Arc Minutes Arc Minutes Arc Minutes 


22 


Tibbs et al. 


0 20 40 


10 20 H/cm 2 


60 80 100 120 140 


12 


14 


K 


16 


18 



-i_i_l_i_i_i_i_l_i_i_i_,_l_,_,_,_,_1 

5 0-5 -10 

Arc Minutes 

Center: R.A. 04:29:56.29 Decl.+54:14:51.7 
10 20 H/cm 2 


0 50 100 150 200 250 300 



5 0-5 -10 

Arc Minutes 

Center: R.A. 04:47:23.41 Decl.+53:03:31.4 



Center: R.A. 04:29:56.29 Decl.+54:14:51.7 
K 


12 14 16 18 



Center: R.A. 04:47:23.41 Decl.+53:03:31.4 


0 50 100 


10 20 H/cm 2 


150 200 250 300 



14 


K 


16 


18 



5 0-5 -10 

Arc Minutes 

Center: R.A. 04:47:57.83 Decl.+53:07:51.2 


o 

< 


5 0-5 -10 

Arc Minutes 

Center: R.A. 04:47:57.83 Decl.+53:07:51.2 



Figure Bl. Continued 


© RAS, MNRAS 000, ??-[I7| 


































